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We study the problem of the crossover from one- to higher-dimensional metals by considering 
an array of Luttinger liquids (one-dimensional chains) coupled by a weak interchain hopping t± . 
We evaluate the exact asymptotic low-energy behavior of the self-energy in the anisotropic infinite- 
dimension limit. This limit extends the dynamical mean field concept to the case of a chain embedded 
in a self-consistent medium. The system flows to a Fermi-liquid fixed point for energies below the 
dimensional crossover temperature, and the anomalous exponent a renormalizes to zero, in the 
00 ' case of equal spin and charge velocities. In particular, the single-particle spectral function shows 

sharp quasiparticle peaks with nonvanishing weight along the whole Fermi surface, in contrast to 
the lowest-order result. Our result is obtained by carrying out a resummation of all diagrams of the 
D ' expansion in t±_ contributing to the anisotropic D ^ oo limit. This is done by solving, in an almost 

5_( , completely analytic way, an asymptotically exact recursive equation for the renormalized vertices, 

"^ ■ within a skeleton expansion. Our outcome shows that perturbation expansions in tx restricted to 

lowest orders are unreliable below the crossover temperature. The extension to finite dimensions is 
C^ , discussed. This work extends our recent Letter [Phys. Rev. Lett. 83, 128 (1999)], and includes all 

mathematical details. 

^ ' PACS numbers : 71.10.Pm, 71.10.Hf, 71.27.4-a 11. 10. Hi 
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I. INTRODUCTION 

J> r According to Fermi-liquid (FL) theoryBB, a quasiparticle is identified by a single dispersive coherent peak in the 

■^ single-particle spectral function describing a particle or a hole close to the Fermi surface (FS). This peak becomes 

'""' ■ sharper when approaching the FS, which reflects the fact that the lifetime of the quasiparticle becomes infinite at the 

^^ ■ FS, while keeping its total weight Z (quasiparticle weight) finite. On the other hand, FL theory fails generically in one 

^^ , dimension, where quasiparticles are not well defined, and the elementary excitations consist of collective charge and 

Qs spin excitations with bosonic properties. In this case, the single-particle spectral function shows two dispersing peaks, 

^s^ corresponding to charge andscin modes. The splitting into two peaks corresponds to the decay of the quasiparticle into 

"Y^ spin and charge excitationtu B, i. e., the spin and the charge of an injected electron move independently with different 

Ci ' velocities. A more important result is the fact that the quasiparticle weight Z vanishes when the FS is approached. 

S . This implies that for k equal to the Fermi momentum kpj where spin and charge energies merge, the spectral function 

I ' does not become a delta function as a function of frequency w, but rather it diverges with a weaker power-law behavior 

' ^ , like w"~^. This refiects onto the behavior of the momentum distribution n{k), which no longer shows a discontinuity 

i^ ' at fc = fci?, but rather a power-law behavior {\n{k) — n(fci?)| ex |fc — fc_F|"). The same exponent appears in the local 

f_\ \ density of states, which vanishes at a; = like w". The exponent a thus characterizes the anomalous behavior of 

• • ■ one-particle correlation functions and it plays the role of the anomalous dimension as in field theory. However, in 

contrast to the usual field-theoretical models (as -theory) , the anomalous behavior of one-dimensional Fermions is 

not universal, since the exponent a depends on the interaction. ID metak kaving these properties take the name of 

$H Luttinger liquids (LL), the name coming from the Luttinger model (LM)LrLl, which plays the role of the "canonical 

model" for 1-D interacting fermions. ■— ■ ■— . 

The interesting question is what happens between one and two dimensional3ilj. Specifically, one can start from a 
D-dimensional array of chains (the interesting cases are, of course, D = 2 or 3), initially uncoupled, and then switch 
on a small tunneling (hopping) amplitude i_L between the chains. The question is when and how does the crossover 
to a normal FL behavior occur? While the question of the crossover from an anomalous LL to a normal FL state 
is a challenging problem per se, there are other reasons why one is interested in this problem. The first two are 
connected to the theory of high-Tc superconductivity. First, it has been suggested that-tJie, normal-state properties 
of high-Tc superconductors may be explained by some kind of two-dimensional LL statellj'EZl. Once a 2D LL state is 
assumed within a_Cu02 plane, it has been suggested that incoherent hoppingj-hetwecn different layers may favor a 
BCS paired stateEa. Secondly, it has become clear from a variety of experimentfO that underdoped high-Tc materials 



are characterized by the presence of charge modulations in the form of one-dimensional stripeala. In these structures, 
the electron dynamics occurs mainly in the direction longitudinal to the stripes, and, t hu s , |it could be effectively 
described by quasi-one dimensional models in which the transverse dynamics is reducedEfcl. The third reason is 
related-to-the existence of several synthetic and natural compounds which can be-considered as quasi-one dimensional 
metalSEl, such as the organic conductors TTF-TCNQ, the Bechgaard saltsN {TMTSF)2X and {TMTTF)2X 
(with X = PFq, CIO4, • • •), or the inorganic chains NbSe^, Kq,^MoO^. A further possibility to study the crossover 
between ID and 2D is to couple a finite number of chains together. The phase diagram of suchJaddjer systems is quite 
rich, and it shows an interesting dependence on whether the number of chains is even or odd o EJ. 

In this paper, we consider the effect of a small tunneling matrix element t±_ coupling the chains. The question is: 
does the electron liquid go over to a FL state for arbitrarily small i_L and sufhciently low temperatures or is there 
a critical value of ij_ below which one has a LL state for arbitrariUtjlew temperatures? This question is related to 
the problem of dimensional coherence addressed by Anderson et alll3't3. These authors suggest that for sufhciently 
strong interaction the system may remain in a LL state for sufficiently small t_L . Clearly, the correct starting point, 
as stressed by these authors, is to consider initially the problem of uncoupled LL and then treat ij_ as a perturbation. 

However, renormalization-group calculations show that t± is a relevant perturbation which means that an arbitrarily 
small ij_ should destroy the ID LL stateEJ. 'Riis can be understood from simple dimensional arguments. Consider 
the LL Green's function t/(x|0) in real spaceS This goes like |x|~^~" at large distances, and thus the Fermi field 
operator ^(x) ex •\/^(x|0) has dimensions [^(x)] = £^jp^^"''^. Therefore, upon integrating over the imaginary time 
r, the perturbation associated with the t±^ term (see Eq. (nl) below), has dimension fy"""'^. This means that each 
term in the perturbation expansion in t±^ carries a term E\"~^ , which diverges at low energies whenever a < 1. These 
divergences signal the fact that the perturbation t±_ is relevant for a < 1. 

Let us consider the energy at which higher-order terms in the tj_ perturbation start to become important (i. e., 
all of the same order). This is given by fy = tj^^''^^"' = te//- This introduces a new energy scale, te//, which 
characterizes, for example, the,(ppssover temperature above which temperature fluctuations cover the effect of t_L and 
the system behaves like a LlJl3'E^. This means that for temperatures T much smaller than Ep but much larger than 
te// the scaling behavior is characterized by the LL anomalous dimension a. For example, the Green's function at 
k = kp scales like w""^ (for w 3> T) in this range. In this temperature region, the system is still effectively one 
dimensional since the effects of t± are washed out by the temperature. Below this crossover temperaturellj and for 
energies smaller than t^ff the effects of t± become important and higher-dimensional coherence sets on. Notice that 
the effect of electron interactions are indeed important in reducing the coherence of the interchain hopping. In fact, 
the crossover temperature is reduced considerably for a > 0, since, in this case teff <?C t_L, and the interchain hopping 
maintains an incoherent behavior down to very low temperatureslij. However, strictly speaking, whether the system is 
a FL, a LL, or something else can be determined in the T — > limit only, since both of them are asymptotic theories, 
i. e. valid in the low-energy limit. Therefore, the important energy region to be studied is fy <C teff- This is the 
nontrivial region, since the behavior here is determined by all terms in the t± expansion. 

For this reason, any perturbative expansion restricted to lowest order is uncontrolled at low energies £y <C ^e//, 
and lowest-order expansions are inconclusive. This is the reason why theoretical results are still contradictory about 
the nature of the ground state in this energy region. Since, as discussed above, this is precisely the relevant region 
for a possible FL behavior, it is worthwhile investigating it in a controlled way. This has been done in Ref. O, by 
considering all diagrams corresponding to the infinite-dimension limit. In this paper, we extend the results of that 
Letter, and provide the details of the calculation 

This paper is organized as follows. In Sec. |^, we introduce the problem of LLs weakly coupled by a single- 
particle hopping tj_ . We discuss the issue of the perturbation expansion in t± , its difficulties, and the lowest-order 
approximations. Next, we discuss the limit considered here, namely the "anisotropic" D —^ 00 limit, and the analysis 
of the asymptotic low-energy regime. Finally, we present an appealing discussion of the analogy of our method with 
the parquet summation and with the renormalization group, and discuss the cases in which the present method is 
controlled. In Sec. |l^, we describe in detail the procedure to carry out the sum of the diagrams leading to the D — > c» 
limit for the self energy Eq. (0) . The idea is to write a recursive equation for the "restricted renormalized cumulants" 
Eq. (^) in terms of the effective hopping T±. In the leading logarithmic order, this gives a set of self-consistent 
recursive equations, Eq. (S), which can be easily solved to a very high degree of accuracy by a power expansion and 
a Fade analysis. In Sec. |^, we discuss the results of this calculation. The most important one is the fact that the 
anomalous exponent scales to zero, i. e., the self energy no longer scales anomalously at low energies. This is seen 
in the spectral function close to the "special" Fermi point c± = 0, which becomes sharper, in contrast to the the 
lowest-order approximation. The quasiparticle weight no longer vanishes at cj^ = in our result. Finally, in Sec. M 
we state our conclusions, and discuss possible extensions of the calculation to the inclusion of spin-charge separation 
and to finite dimensions. 

Due to the absolute novelty of our procedure, we considered that the reader would benefit from an inclusion of all 



details of the calculations, so that any one could follow and repeat our steps without difhculties, and possibly extend 
them to some other cases. The calculation is transparent, as it is almost completely analytic except for the Pade 
solution of the recursive equation described in Sec. |^. In order not to burden the bulk of the paper, we deferred most 
of these calculational details to the appendices. 

II. THE PROBLEM: FROM ONE TO HIGHER DIMENSIONS 

We consider a. D' {= D — l)-dimensional hypercubic array of parallel one-dimensional chains (i. e., the total 
dimension is D). We consider here the case of equal spin and charge velocities, since it allows for crucial simplifications 
in the calculation. Since we are interested in the effects and in the fate of the anomalous exponent a, we believe that 
spin-charge separation should not play an important role. The chains are labeled by the D — 1-dimensional coordinate 
x± along the hyperplane perpendicular to them, while the coordinate along the chains is called xu . The Hamiltonian 
we want to study has the following formS: 

H = '^Hll{x±) + ^ t±{x±-x±) '^ dx\\ 'iIjI^^{x\\,x±) ipr,cr{x\\,x±) , (1) 

x± x±x±' rcr 

where tpr,a{x\\,x±) [tpl ^{x^\,x±)] is the destruction [creation] operator for a right- (r = +1) or left-moving (r = —1) 
fermion at the position xh along the chain x± with spin a. Moreover, Hll{x±) is the Hamiltonian for an (uncoupled) 
LL in the chain x±. Since we are interested in Jcwfc-energy properties we can just take for Hll{x±) a Luttinger 
model, characterized by its parameters a and w^LrQ'B (since we neglect spin-charge separation), which will depend 
in a nontrivial way on the bare parameters of the microscopic chain Hamiltonian. However, we are not interested 
in this dependence here, and we just take these parameters as our starting point. In Eq. (m, t±{x± — x±') is the 
amplitude for the hopping of an electron from chain x±' to chain x±, where, as usual, we have assumed that neither 
the Xji coordinate, nor the direction r are changed by the hopping. Moreover, one can restrict to the case of an 
hopping between nearest-neighbor chains only. Inclusion of an hopping with finite extension in the xn direction, or 
of a next-nearest-neighbor hopping in the x± direction is straightforward. However, it is not expected to change the 
low-energy results. With tj_ = 0, the problem can bo-solved exactly, as the ground state is given by the product 
ground states of the LM in each chain, which is knownQ 13. 

Knowing the exact solution of the t± = problem, one can envisage carrying out a perturbative expansion in 
powers of t± , as t± is small. This is, however, not without complications, as Wick's theorem does not hold for 
the t± = ground state, since the LM, although exactly solvable, contains electron-electron interactions. A similar 
problem occurs for the expansion about the atomic limit of the Hubbard model, whereby one first solves the single- 
site problem exactly and then expands in powers of the hopping t. A diagrammatic formulation for this problem 
was introduced by Metzner in Rcf. |4[ and further discussed in Ref. |5|. It consists in carrying out a linked-cluster 
expansion, where an arbitrary (even) number p£ lines (2n) can join into one dot. This dot is associated with the exact 
n-particle cumulant of the single-site problemc3. 

This method has been extended to the problem of expanding about the LLs in Ref. |l^. The diagrams contributing to 
the expansion are the same, the only difference being that each line is now labeled by the extra variable xy (intrachain 
coordinate) and r (for left- or right-moving fermions), besides spin a, and imaginary time t. Actually, this method 
turns out to be more appropriate for the present problem rather than for the Hubbard model. Indeed, in the Hubbard 
model, one expands about an highly degenerate i = ground state, which is not the case in our problem of coupled 
LLs at t_|_ = 0. Alternatively, one can use the diagrammatic rules in momentum space, for which each line carries 
an intrachain momentum fcii, a Matsubara frequency uj, and an interchain momentum kj_, as well as indices a and 
r. Besides this modification, rule 2 of Ref. p3 for calculating the Green's function remains the same. A set of these 
curious diagrams, contributing to the Green's function, are shown in Fig. |l|. The building blocks of the diagrammatic 
expansion are (i) hopping lines connecting nearest-neighbor chains (say a;_L2,a;_L2) associated with t±{xj_i —xj_2)j ^^id 
(ii) "dots" with n entering and n leaving legs, associated with the n-particle cumulant of the single chain. The latter 
can be readily evaluated, at least for low energies, since one knows the exact solution of the Luttinger model and of 
its correlation functions (cf. Sec. ^). _. 

Boies et al. used a functional-integral method to obtain an expansion in t± about the LLIIj. Although their 
formulation allows, in principle, for an expansion to any order in tj_ , in practice one can just get the first few 
orders. Our method provides a systematic diagrammatic formulation of this expansion to any order. The advantage 
of a diagrammatic formulation is that one can choose a class of diagrams to sum over, according to some physical 
guidance, without being restricted to the few lowest-order terms. This is particularly important for the model at 
study, since, as discussed in the Introduction, each power of t± in the perturbation carries a term £\\"~ , which 
diverges precisely in the important region. Thus, one cannot reliably restrict to a finite number of diagrams. 
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FIG. 1. Diagrammatic expansion in t± of the single-particle Green's function Q (grew box). A directed line connecting two 
chains x±i and x±2 gives a contribution t±{x±i —x±2)y or t±{k±) in momentum spaca23. A dot with n entering and n leaving 
lines contributes a factor t/° (n-particle cumulant of the uncoupled LL, see Sec. En]), (a) Example of single-particle irreducible 
and reducible contributions to Q. (b) Dyson's equation for Q in terms of the inverse-self-energy F (gray disk), (c) Example of 
diagrams contributing to F. (d) Self-consistent diagrams contributing to F in the D ^ oo limit. The self consistency is due to 
the presence of the full Q in the internal lines of the loop. 

Some diagrams contributing to the expansion of the Green's function Q (gray box) are shown in Fig. |l|a. As in 
conventional perturbation theory, one can consider the function T obtained by the sum of irreducible diagrams, i. e., 
the ones which cannot be separated by cutting a single line (see Fig. |l|c) . One then obtains a Dyson- like equation for 
Q as a, function of F (Fig. [Ip) of the forn£3 

C?(k) = [F(k)-i-i^(fc^)]-' . (2) 

Notice that F"-'^, and not F, appears in the inverse Green's function contrary to standard perturbation theory. For 
this reason, we call F inverse self energy. 

The lowest-order approximation for F (the "dot": a in Fig. |l|) copesponds to taking T — Q^, the Green's function 
of the isolated LL. This gives for the total Green's function Eq. (||)E3 

g{k) = [g%k)-'-t^c^]~\ (3) 

This expression is a generalization of the Hubbard I approximation for the case of an expansion_about the LL. Eq. (g) 
has been first obtained by Wen via a different procedureEZi, and reobtained by Boies et al.ll^ within a functional- 
integral method. This approximation, which we will refer to as "LO ", is also called "single-dot", "RPA", "Wen's", or 
"Hubbard I" in other papers. For a < 1, the effect of the interchain kinetic energy t±c± is to change the branch-cut 
singularity into a true quasiparticle pole (cf. Ref. Bq) for all k points close to the FS, except for those kj_ points for 
which c± = (for example, for D = 2 these are k± ~ ±n/2). In particular, the positions of the poles ioi lu — 
identify the new FS, which acquires a dispersion of the form /cy p{k±) ex {t±c±)^'^^~°'\ i. e. it is rediiced with respect 
to the noninteracting case, where one would have fcy p,{kj_) ex tj_cj_, but not completely suppressedll3. For the sake of 
completeness, we discuss the main results of this approximation in Sec. ^. 

Since the branch cut are shifted into poles, this approximation gives a FL along the whole FS except close to the 
c± — region. This can be also seen from the quasiparticle weight Z, plotted in Figjj (dashed line), which vanishes 
for c± = 0. For this reason, the quasiparticle peak is quite broad in this region, as can be seen from Fig.H. However, 
as discussed above, this result, being restricted to lowest order is uncontrolled in the region fy <^ teff and one should 
sum an infinite series of diagrams in order to get reliable results. Since it is not possible to sum all diagrams in 
the expansion, we want to select a workable subset of diagrams according to some physical limit in order to avoid 
an arbitrary choice. Specifically, we consider the series given by the diagrams indicated in Fig. mi, corresponding to 
the large-dimensioa limit (D — > oo). The D — > oo procedure adopted here is different from the standard dynamical 
mean-field theory LJ, since our system is strongly anisotropic, as the hopping in one (in the ||) direction is not rescaled 



by the usual 1/^/^ factor and is much larger than in the other D — 1 (_L) directiong23. In analogy to the standard 
D —> oo methodo, where one has a single impurity embedded in a self-consistent medium, our I? — > oo system 
represents a 1-D chain embedded in an effective self-consistent medium. As a|-consequence, the self-energy is local 
with respect to the _L coordinates but has a nontrivial dependence on the || onesEj. We believe that this is the correct 
starting point to study the crossover problem, since, in this way, one treats the one-dimensional problem exactly and 
includes the coupling to the other chains by an effective dynamical mean field. 

Even summing all the D = oo diagrams is an impossible task. Nevertheless, since we are interested in low-energy 
properties, we can restrict to the leading singularities in each diagram. It turns out convenient to rewrite the power 
expansion in terms of the dressed hopping Tl (indicated by a dashed line in Fig. m). This is very similar to the 
skeleton expansion in conventional perturbation theory, where self-energy insertions are removed. The advantage 



is that the scaling behavior of the effective hopping (cf. Eq. (C12)) exactly cancels the power-law divergences of 
the diagrams, and each term of the perturbation acquires the same scaling as a function of the energy, and only 
logarithmic divergences are left. 

The procedure of summing just the leading logarithmic divergences-is similar in spirit to the sum of the leading 
divergences in the parquet series, which was introduced by the Russiancd and by the Frcnchc3 school in order to study 
the instabilities of various one- and higher^dimensional electron systems. This method is equivalent to the one-loop 
renormalization-group (g-ology) approaches, and it actually gives a rigorous background, as well as a systematic 
formulation for the extension of the g-ology method to higher dimensions. In our case, this corresponds to considering 
the quantity I = a\og{t^ff/£i\) to be of order one, mid thus taking all orders in I, while considering a small. 

Similarly, in the parquet summation, or g-ologycS, the small parameter is the bare interaction vertex go and one 
sums all powers of go log -^^ , uj being the characteristic energy scale. The sum of this series gives the renormalized 
interaction vertex g{uj/Ep) which thus acquires an energy dependence. Within the renormalization-group picture, 
the energy-dependent interaction vertex is interpreted as an effective interaction acting on an effective low-energy 
subspace, i. e., on a subspace in which high-energy modes are integrated out. Whenever the interaction vertex 
scales to zero, this signals that the effective low-energy theory describes non-interacting electrons, i. e., the theory 
is asymptotically (infrared) free. As a consequence, the exponents of correlations functions are mean-field like and, 
in the case of fermions, the system is a Fermi liquid. On the other hand, when a vertex diverges, no controlled 
prediction can be made about the low-energy behavior of the system, since the perturbative approach breaks down 
for sufficiently low energies, even when go is small. In this case, the divergent vertex signals an instability towards 
some kind of broken-symmetry state. 

In our case, the role of the interaction vertex is played by the anomalous exponent a. The bare a is the correlation 
exponent of the uncoupled set of Luttinger liquids. Switching on the interliquid hopping t± produces a renormalization 
of the exponent. This renormalized exponent is obtained by looking at the low-energy behavior of the self-energy in 
the coupled-chains system. Similarly to the f;-ology case, our result, obtained by summing the leading logarithmic 
divergences, is thus controlled if (i) the starting (bare) value of a is not too large and (ii) a scales to zero for low 
energies. The first (i) requirement is easy to fulfill, since for most interesting systems a is quite small For example, 
for the Hubbard model a < ^, where the equal sign holds for an infinite value of the on-site interaction U. Larger 
values of a are obtained by increasing the range of the interactiorO. This is another reason why our approach is 
more convenient than a weak-coupling expansion in U: while our calculation makes sense also for very large (bare) U, 
for which a is still small, the weak-coupling renormalization group is not justified for U larger than the bandwidth. 
An estimate of the maximum value ac of a, for which our calculation is justified is given in Sec. M. The second 
(ii) requirement can be only checked a posteriori. The main result of this paper is that indeed point (ii) turns out 
to be satisfied, as a scales to zero for energies smaller than t^ff . Thus, our procedure of restricting to the leading 
logarithmic divergences is controlled, unless one starts from a model with a too large value of a. 

III. ANISOTROPIC D ^ oo METHOD 

In this section, we carry out the piim of the D —^ oo diagrams for the inverse-self-energy. In the D —^ oo limit, the 
inverse-self-energy r(xo) is _L-localEl, and is obtained as the sum of the loop diagrams in Fig. ||a (equivalent to the 
ones of Fig. |l|d) as 

r(xo) = e°(xo|0) + Y, — / [n '^^y^ ^'^'^' ^i(-Xfc>0)J g°iyo + xo, • • • ,y™ +x,„|yo, • • • ,y™) , 

m=l "^' •' fc=l 
oo -. m 

= g°(xo|0) + ^(-l)™ / \Y[d^ykd^Xkm-^k,0)\g°{yo+^o,---,ym+^,n\yo,---,ym), (4) 

m=l "'Um k=l 
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FIG. 2. (a) Diagrams contributing to the inverse-self-energy F in the D = oo hmit within an expansion in the dressed 
hopping T± (dashed hne). (b) Dressed hopping and its diagrammatic expression in terms of the bare hopping t± (fuU hne) and 
the Green's function. Other conventions are as in Fig. hi. 



where in the last line wc have exploited the symmetry for exchange of the coordinates I, ■ ■ ■ ,m and restricted the 
integration to the region |xi| > |x2| > • • • > |xm| indicated by "1 | m". The corresponding factor ml is then canceled 
by the symmetry factor 1/m! of the diagram. In Eq. (Q), GciVo ' ' ' ymlyo ' ■ ' ^m) is the m + 1-particle cumulant of 
the uncou pled LL, i. e. the c onnected part of the m + 1-particle Green's function Gciy'o ' ' ' y'm\yo ■ ■ ■ ^m) defined in 
Eq. ( All| ) (see also Eq. (D23) for the definition of cumulants in terms of Green's functions). In particular, for m = 
the single-particle cumulant Q^ coincides with the Green's function G°, as there are no disconnected parts. Moreover, 
Tl(x, a;_L = 0) is the dressed hopping written in real space, which is calculated in Sec. O. 

We are interested in the dominant low-energy behavior (fy <C ie// corresponding to |xo|ie// 3> 1) of correlation 
functions and thus we can restrict to the leading logarithmic divergences in the loop integrals ( Eq. (^)), as discussed 
in Sec. |l|. Let us estimate this leading contribution. If, as a first step, one neglect s the self-consistency of the Green's 
function and dresses the hopping T± with the bare Green's function only ( Eq. (C12)), one can see that the leading 
contribution of a 77i-loop term in Eq. (^) has the form tj^(xo|0) x (alog Ixolte//)"^™. Indeed, one "alog" term arises 
from each integration of the "center of mass" coordinates y^ (cf. Sec. 0), another "a" from each T^, due to its 
real-S|pace structure (cf. Eq. (C12)), and a "log" comes out for each integration of the "relative" coordinates ( 
Eq. (I). 

Even summing up "just" the leading logarithmic divergences of the integrals in Eq. (^) is a tough task. To do this 
we proceed in several steps. First, consider that some integration regions in Eq. (||) can be left out, as they don't 
contribute to the leading logarithmic divergences. Specifically, in addition to the region_lxi| > |x2| > ••• > \xm\ 
(called 1 I m), to which we restrict by symmetry, we can further restrict to the regionEj |xo| > |xi|, and |xp| < 
min{\yq — y,|, |y^ — yr|, \y'q — y'r\) for each p > q,r (of course, q ^ ?■), and y^ is defined as y, -|- x^. The fact that 
the leading logarithmic contributions only come from this integration region, which we will call "0 -IJ- ttt." , is proven in 
Sec. I 



For convenience, we introduce the "restricted renormalized cumulants" (RRC)Cj, defined only in the region "0 JJ- 771" 



as 



Gc{yo + xo, • • • , y™ + x™|yo, • • • , y„0 
= Gc{yo+^o,---,yrn +x™|yo,---,y,„) 

- / (fx,n+l (fym+1 Tl(-X„i+i,0) 

-'OJLm+l 



(5) 



'xQciyo +xo,---,y„i+i +x„,+i|yo,---,y„i+i) . 

Comparing Eq. (^ and Eq. (|^), it is straightforward to verify that r(xo). is given by the single-particle RRC, 
ec(xo|0). 

We thus proceed by evaluating the integrals in Eq. (M). An important point, which we will show below, is that, at 
the leading logarithmic order, the m + 1-particle cumulant is renormalized by a multiplicative factor, which depends 
on the absolute values of the relative coordinates |xi| only. More precisely, the RRC can be written as 



Qc{ya + xo, • • • , ym + x„|yo, • • • , y,„) 

= Tm{lo, ■■■ ,lm) ^c (yo + Xo, • • • , ym + X„i|yo, • • • , Ym) 



(6) 



where the J^m is the renormalization factor, which we have written in terms of the logarithmic variables It — 
a\og{\xi\teff)- We can thus first carry out the integration over the "center-of-mass" coordinate y,u+i in Eq. (||) by 



simply considering the effect on the bare cumulant, as the renormalization factor does not depend on ym+i- This 
integral is quite involved, bui-|its leading logarithmic contribution can be calculated analytically. This is carried out 
in Sec. ^ where one obtains^ 

/ rf^ym+i ^c(yo +xo,---,ym+i +x„,.+i|yo,---,y„i+i) (7) 

OiJ.m+1 

= 27r t;°(yo+xo,---,ym +x„|yo,---,ym) 

m . ^ 

X|x„+ip ^°(x„j+i|0) ^a (log|Xj| -l0g|x,„+i|) (^r,,r,„ + i + 7; (^r,-, 



J = l 



-»-™+l 



After carrying out the integral over ym+i we carry out the integration over the "relative" coordinate x^+i, which 
includes a sum over r and ctH. Inserting the form Eq. (g) and the result Eq. (0) into Eq. (||), and dividing both 
sides of the equation by C/0(yo + xq, • • • , y™ + x^lyo, • • • , y™), one obtains 

^m(?o,---,^m) = l-2^ f d2x„+iri(-x„+i,0) |x™+i|2 5,"(x„+i|0) (8) 

ie//~^<|x™ + i|<|x,„| 



s 



-I-m + l I J 



where the lower limit of integration for |xm+i| is due to the fact that Tl changes its behavior in the region Ep < 
|xm+i| < teff"^ (Ref. [47| ) , and, thus, there is no logarithmic contribution here, and the upper one i s due to the 
restriction IJ- to + 1 in Eq. (|5|). Inserting the asymptotic expression for the dressed hopping Eq. ( C21 ) in Eq. (|g). 



one can carry out the integration over x„j+i in circular coordinates, and obtain the recursive self-consistent equation 
for T^ 

^m(?0, • • • , ?m) = 1 + 2(1 + 5) / dlm+l V {Ij " Ira+l) (9) 

X ^m+l{lo, • • • , Im+l) [^o{lm+l) + ^o{lm+l)\ ■ 

From Eq. (m, it is obvious that Jv, depends on just two variables, namely I = Ig + ■ ■ ■ + Im-i, and Z„j. With this 
redefinition, and renaming the integration variable Im+i to /', Eq. ffl) can be reduced to 



TmilJrn) = 1 + 2(1 + 3) / dl' [I + kn - [m + I) I'] Tm+lil + l^ J') [MH + ^oiH] ■ (W) 

JO 

Eq. (no) is a self-consistent equation, since Tq — l/Trn=o [ Eq- (ClTp], which depends on the Tm, to insert on the 






r.h.s.. We have not been able to find an analytic solution to Eq. (|10|). However, by expanding in powers of of the 
variables h one can write a recursive equation for the coefficients of the expansion of the functions J-'m up to a rather 
high order with a moderate numerical effort. This procedure is described in detail in Sec. R 

We have evaluated the coefficients of ^o up to the 42*^ order in L A Mathematica program has allowed us to evaluate 
these coefficients in a rational form, which is particularly recommended for a Fade analysis. A straight summation 
of the series is not recommended, since its convergence radius seems to be rather small (of the order unity), while 
we need the asymptotic behavior for large I. Nevertheless I — alog(|x|ie//) is restricted to the neighborhood of the 
real positive axis, and a Fade analysis shows that the poles are either away on the complex plane or on the negative 
real axis. A Fade analysis is thus the most appropriate procedure in order to determine the large-Z behavior of the 
function J-o{l)-, which also gives the asymptotic behavior of the inverse self-energy r(a;). The results will be presented 



and discussed in Sec. W 



IV. RESULTS AND DISCUSSION 



As shown in Sec. ^ the solution of Eq. I fivi ) gives J-o(Z) ~ e'^' for large /, where the exponent c turns out to be 
essentially equ al to 1 (within about 10"'' of accuracy) in both cases with and without spin. Introducing this result 
and Eq. ( |A12 ) in the expression for the inverse-self-energy ( Eq. (g) with m ~ 0) yields 




FIG. 3. Fermi-surface dispersion fcy as a function of the off-chain kinetic energy c±_ (in units of ix ) for the coupled 
spinful Luttinger hquids [ Eq. (hi)] with bare LL exponent a = 1/4. Our D ^ oo result (solid line) is compared with the LO 
approximation Eq. (g) (dashed). 



r(x) = (?0(x|0) ^o(alog(|x|V/)) ^ gO(x|0) (|x| t^j^T « W/"/|x| 



(11) 



i. e. the anomalous exponent a exactly cancels out in the asymptotic behavior of F!. The same thing happens in 
momentum space. From Eq. (C16) one notices that the (asymptotic behavior of the) renormalization function is the 
same in momentum space, provided one replaces |x| with l/|k|. Thus, for low energies we obtain for the right- moving 
component {r = -\-l) 



r(k) ^ gO(k) te//" |kr" « t 



eff 



[iLO ■ 



fl2) 



where we have used Eq. (C8). 

The Green's function of the coupled system is given by the Dyson equation Eq. (0). Taking the result Eq. (^3), 
one can readily notice that the Green's function now has poles at ito — k\\ oc tef f°'t±{k±)^ i. e. even for c±_ = 0, in 
contrast to the LO result, where a branch cut was present. In particular, at the FS [ilo ^ + iO^) and for c±^ = 0, 
our result becomes asymptotically exact, as |k| vanishes at the pole. 

Let us look at the FS more precisely. This is the curve kn „{c±_) parametrized by the Fermi momentum as a function 
of the _L momenta, and is determined by the solution of the equation r(fcj| p(cj_), ito — Q-\- iO^)^^ = t^c^. Obviously, 
Eq. (12) gives fcy ^(cj, = 0) = 0. In Fig. g we plot the FS curve for other values of cj^, and a = 1/4 in the case of 
particles with spin. We compare our result (full line) with the LO result (dashed line). For small c^_ , our result 
gives a regular behavior /cm p{c±_) oc t^ff c± in contrast to the lowest-order result, which gives a flattening of the FS 

at cj^ = 0, due to the behavior fc||^(c_L) oc tg// c±^^^^~°'\ 

The quasiparticle weight Z{c±) at the FS is given by the inverse of the coefficient of the linear term in lo in the 
inverse Green's function, more precisely, Z{cj_)~'^ = ^l/5(fc||p(c_L), ia;)jj^^o-i-io+- We have plotted Z as a. function 
of c_L for the case with spin in Fig. 0, again compared with the LO approximation. Moreover, in order to show 
the importance of summing the infinite series of diagrams, we have included the result obtained by truncating the 
D ^ oo series (Fig. ga) at the first loop, by still taking the self-consistently dressed hopping as internal line. For 
small c± , the lowest-order result (dashed line) gives a Z vanishing as Z{c±) oc (ij^cj.)"/'"'^""^ thus yielding poorly 
defined quasiparticles around c± = 0. Inclusion of the first loop (dotted line) gives a vanishing Z too. Therefore, self 
consistency is not enough to restore the FL behavior. Our result, instead, yields a finite Z for c± — *■ 0, as can be seen 
from the figure (solid line). The correct FL behavior is thus recovered on the whole £S, including the regions c± — 0. 

These results can be more concretely seen in the spectral function for small cj^ES. This is plotted in Fig. for 
different cj_, and for fcy = 0. The figure shows a well-defined dispersive quasiparticle peak, which becomes sharper 
by approaching the FS, as should be the case for a FL. The dispersion as a function of c± is a clear indication of 
higher-dimensional coherence. For comparison, in Fig. 0, we have shown the LO result. As one can see, the peak is 
dispersive too, but much broader and lower (notice the different scale). Moreover, a closer inspection shows that the 
quasiparticle weight decreases by approaching the FS, which is consistent with what we have shown in Fig. ^. 




FIG. 4. Quasiparticle weight Z as a function of the off-chain kinetic energy c± with the same conventions as in Fig. p^ . In 
addition, we show the result (dotted line) obtained by partially improving on the LO approximation, i. e., by including the 
first self-consistent loop for the inverse self-energy of Fig. nld. 
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FIG. 5. Spectral function A{lo) of the coupled spinful Luttinger liquids for difi'erent values of Cx, 



0, and a = 1/4 



from our D ^ oo result. For the sake of clarity, the different curves are shifted vertically by steps of 10. They correspond to 
c± — 0.05,0.1,0.15,0.2,0.25,0.3,0.35, from bottom to top. Notice the sharpening of the peaks upon approaching the FS at 

c± ^0. 
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FIG. 6. Spectral function A(lS) from the LO approximation for the same parameters as in Fig. pi except that the curves 
are shifted by 1. The peaks sharpen upon approaching the FS, but the quasiparticle weight vanishes. 
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FIG. 7. Spectral function of the isolated Luttinger model (with equal charge and spin velocities) for a = 1/4 and k\\ — 



We want to study the spectral function even for fcy ^ 0. To understand what happens, let us first look at the 
spectral function for the LMau (without spin-charge separation), which we plot in Fig. ^ for fc|| = 0.2. From the 
figure, one can readily recognize the two nonanalicities at w ± fcy . For lu \ +A:|| one has in fact a divergence like 
{u; — fcj|)"/^^^, while for lu /" — fcy the spectral function vanishes as (/cy + w)"/^. The power-law divergence instead 
of a pole at w = fey is due to the fact that the point, w = fey, where the inverse Green's function, 1/5", of the LL 
vanishes, is not a simple zero but a branch cut. Between ztfey the spectral function of the LM is identically zero, as 
the Green's function has neither cuts nor poles here. At fe|| = the two nonanalicities merge in a single power-law 
divergence uj'^~^. 

Within the LO approximation, Eq. (p[), the zero of \/Q is shifted away from the branch cut. Thus, an isolated 
quasiparticle pole appears in the region —fey < w < fey on the real axis. This pole is always present for any cj^ 7^ 
(see Sec. pi). The pole removes spectral weight from the peak at a; = fey which is now no longer a divergence. 

In our D ^ 00 result, the situation is similar. However, Fig. |9| shows that in this case the two singularities at ±fcy 
loose much more spectral weight in favor of the pole. This is another reason why the quasiparticle weight remains 
larger within our result, as shown in Figs. H, and 0. For c± — Eq. (0) docs not have quasiparticle poles, while our 
result yields a pole with nonvanishing weight at the FS even in this case. The reason for that is due to the different 
behavior of the quasiparticle weight, as shown in Fig. 0, and by the fact that the scattering rate does not vanish fast 
enough for Eq. (g), while it vanishes faster than linearly within our result, as discussed in Ref. ^ (cf. Fig. 2c of that 
reference) . 
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FIG. 8. Spectral function within the LO approximation for the coupled spinful Luttinger liquids with a = 1/4, fc|| = 0.2, 
and c_L = —0.2. In order to make the quasiparticle delta function visible, we have added a small imaginary part ~ 3.0 10""". 
Due to the proximity of the singularity, the peak actually becomes broader. 
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FIG. 9. Spectral function for the D ^ oo result with the same parameters as in Fig. |8[ Notice the much larger transfer of 
spectral weight from the singularities to the quasiparticle pole. 
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V. CONCLUSIONS 

In conclusion, we have studied the problem of the crossover from one to higher dimensions for fermionic systems, 
when Luttinger liquids are coupled by a small hopping t± . Specifically, we have concentrated on the region below the 
single-particle crossover temperature, £n <C ^e//, which is the one relevant for the dimensional crossover. We have 
carried out an expansion in powers of tj_ , and summed the self-consistent series of diagrams (Fig. ITId) corresponding 
to the anisotropic Z? — > oo limit. Our result shows that the LL exponent a renormalizes to zero for energies smaller 
than the single-particle crossover temperature te//- The system thus flows to a FL fixed point with mean- field like 
exponents. This is seen, for example, in the self-energy, which now scales linearly as a function of frequency and 
momentum, in contrast to the LO approximation Eq. (f^, where the self-energy still scales anomalously like |k|"'^~" . 
As a consequence, well defined quasiparticles are recovered along and in the neighborhood of the whole FS, in contrast 
to the result Eq. (j^), where the spectrum is incoherent for small c± . 

We have shown the importance of including an infinite series of diagrams, in order to give reliable results in the 
region fy <C ^e//- Even introducing the first loop of the diagrammatic expansion (in Fig. |l|d) does not give the correct 
result, as shown in Fig. 0. This shows that not even a self-consistent calculation is sufficient. This is the reason why 
previous theoretical results, restricted to lowest orders, are still contradictory about the nature of the ground state in 
this energy region. 

These results have been obtained for the case of equal spin and charge velocities. In fact, we believe that the 
scaling behavior of the anomalous exponent a found here is universal and should not be affected by the inclusion 
of spin-charge separation. Nevertheless, an extension of the present calculation to the case of LLs with different 
velocities could be interesting, first, in order to check this fact, and second, in order to verify whether spin-charge 
separation scales as well to zero in higher dimension like a, or not. 

The imaginary part of the self energy, ~Im F^^, needed to evaluate the, spectral functions in Figs. || and || has 
been determined by analytic continuation of the asymptotic form Eq. (n2[)c3. However, one should mention that our 
calculation, restricted to the leading divergences, yields reliable results for ImT~^ at small values of uj and fcy only. 
On the FS and for large c± , A:|| is large too. Thus, for large c± , we cannot state with certainty whether corrections 
to ImT~^ beyond the leading divergencies vanish fast enough upon approaching the FS or not. Arguments similar to 
the one of ordinary perturbation theorycj cannot be extended to the present case, due to the momentum dependence 
of the vertices in the tj_ expansion. A hint can be possibly obtained by explicitly evaluating numerically the first few 
loops in Fig. nld, without restricting to the leading divergences. 

In principle, we cannot say whether our result is valid also for the physical cases of finite dimensions, and, in 



particular for D = 2 ot D — 3. However, as we have shown in Sec. C 2, the non-_L-local dressed hopping Tx(x, x± ^ 0) 
vanishes faster than the _L-local one Tl(x, 0) for large |x|. Non _L-local contributions are thus irrelevant and one may 
try to extend the present result to finite dimensions. However, there are still J--local diagrams of order 1/D, (for 
example, if one takes the diagram 7 of Fig. m and replaces all internal line with a local T±), which may spoil this 
result. It might be interesting to consider an expansion about the present D = 00 result, and consider the irrelevance 
or relevance of such diagrams, and give predictions about a possible critical dimension Dc, above which the results of 
this paper hold. For example, this could be doiie,in order to study the critical behavior in the neighborhood of the 
transition to the two-particle regime at a = Q;2pli3, where Ugp ~ 0.41 (0.62) for spinless (spinful) electrons. 

In Sec. ||, we have already discussed that our "renormalization-group-like" result holds for a smaller than a certain 
Uc- Although we cannot determine ac exactly within our approach, we can estimate it, e. g., by the value of a for 
which the spectral function becomes negative in some regions. This criterion gives ac ~ 0.50 for the spinless and 
ac ~ 0.33 for the spinful case. 

Another question is the contribution of the shifted poles |k|, c± ^ 0, which turn out to be irrelevant in the present 
case (cf. Sec. O). However, these poles may give important contributions in lower dimensions. Indeed, these poles are 
the one giving rise, in some conditions, to the well known nesting or superconducting instabilities at selected regions 
of the FS. 

The author thanks W. Hanke for useful discussions. Partial support by the BMBF (05SB8WWA1) is acknowledged. 

APPENDIX A: MANY-PARTICLE CORRELATION FUNCTIONS OF THE LUTTINGER LIQUID 

For the sake of completeness, and in order to fix our notation, we give here the expressions for the n-particle 
Green's functions of the LM in real space. To our knowledge, their explicit expression, although known, has not been 



reported anywhere else. In Sec. A 1, we discuss the scaling behavior of the diagrams in the t± expansion. 



A generic _L-local n-particle Green's function is defined aacn 

^'^-••'^-(Xi, • • • ,X2„) =< Tri^^'i^l) ■ ■■i/'"{^2n) > , (Al) 

12 



where Tt is the imaginary-time ordering operator, ip'^iyi) destroys (for d = —1) or creates (for d = +1) a fermion at 
the point x (which includes r and a). In order to extract the t± = cumulants Q^ of the isolated LM, to be used in 
Eq. (Q), we first need the (disconnected) Green's functions Q^ . These can be written as 



iQ,di---d2'n 



(xi, • • • , X2„) = (27ra) "r/ri ■■■Vr 



(A2) 



2n>ii>J2>0 



n [p..^..jx,,-x,j]^''-'- 

;ii>J2>0 

on constraints X]i=i '^i ~ 
fulfilled, otherwise Q^ = 0. Here, a is a short-distance cutoff (a c>c ■gM.__The Klein factors "q^ obey anticommutation 



rules {rfri '7r'} = ^ry and account for the fermionic anticommutationsEjH. From now on, we will set a to unity, unless 
otherwise specified. 



This holds whenever the particle- and momentum-conservation constraints X]i=i '^i ~ 0, and X)i=i '^^''i ~ are 



The functions P in Eq. (A2) can be written as 



Pr.roA^) =7^(x)-9('^l r2,K,) ^i 

where the exponent q{ri, r2, Kp) is given byt3 



q{r,Kp) = 



\i^Kp + ^) 



(f-A(x)) 



for 5' = 1 



i [r{Kp + 1) + -^ + 1] for 5 - 2 



I + a for r = 1 
B for r = — 1 



where the LL exponent a is related to Kp via 



Kp + l/Kp - 2 
2 S* 



and 



Here, 



and 



X\\ = -g-, T = — , 



and 



B = 



_ II Kp - Kp 



25 



7?.(x) = — i/cosh^ a;!! — cos^ t , 



yl(x) = arg(tanha:|| +itanr) , 
At zero temperature T = 4 = 0, Eq. (A7) and Eq. (A8) become 



^(x) 



^'^(x) iarg(a;||+iT) _ 



a 



Xll -H JT X • V 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



where we have introduced the complex vector v = (I,z), allowing for a compact expression. These expressions are 
valid for |x| ^ a and need a short-distance cutoff for |x| ~ a. The cutoff prescription for the LM amounts to replacing 
X|j^ -I- T^ with a;||-^ + (|r| -|- a)^. However, it turns out convenient to adopt a "rotation symmetric" cutoff obtained by 
replacing X||^ -I- r^ with xy'^ + t^ + a^, or by setting 7?,(x) = 1 for |x| < a. The low-energy results, obviously, don't 
depend on the specific choice of the short-distance cutoff. The advantage of setting equal spin and charge velocities is 
clear at this point. Without this assumption, the correlation functions would not be invariant under rotation in the 
(x||, r) plane, which would have made the calculations more difficult. 
In conformity with Ref. |Tl|, we define 
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a"(y[„---,y;„|yor 






(All) 



and Gciy'oi ■ ■ ■ lymlyo: ■ ■ ■ TYm) as the corresponding cumulant (or connected Green's function) to be inserted in the 
diagrammatic expression Eq. ( |4|). 

As an example, we use Eq. ( |A2| ) to evaluate the single- and the two-particle Green's functions (here, we indicate 
explicitly the indices r as 1 for r = +1 or 1 for r = —1). The single-particle Green's function reads 



a°(xl|01) 



27ra 



, |-l-a g-»argx.v 



while the two-particle Green's function for right-moving particles reads 

^0. 1 , Ti M M, g°(yii|yii) g"(y2i|y^i) g"(yii|y^i) g"(/ii|y2i) 

y (yii,y2i yii,y2i) = j^, — r, — ^. ^„, , .. ,.. ■ 

^"(yii|y2i) ^"(yii|y2i) 

On the other hand, the two-particle Green's function for mixed right- and left-moving particles reads 



a"(yii,y2i|y'ii,y^I) = a"(yii|y'ii) a°(y2i|y^i) 



lyi-y-jl |yi ~y2i 
|yi -y2| |yi -y2l 



-B 



(A12) 



(A13) 



(A14) 



1. Scaling behavior of diagrams 



From Eqs. ( Al 1 , A2 , A3 , A9 , A5 ) one can easily extract the scaling behavior of Green's functions for a homogeneous 
rescaling of the coordinates X; -^ Ax^ . 



gO(Ay;„ • • • , Ay;„i|Ayo, • • • , Ay„_i) = A-^^+i) (i+")e°(y^, • • • ,y;.i|yo, • • • ,y„-i) , 



(A15) 



i. e., an n-particle Green's function (and a cumulant too) scales like n one-particle Green's functions in real space. 
Going back to the diagrammatic formalism, Eq. ( A_15) gives the scaling behavior of a vertex with 2n legs. In addition, 
each internal line, associated with a tj_ term, contributes an integration over r and xy, i. e. a factor A^. Let us now 
consider a order- iV diagram {N internal lines), with E external lines. Each internal line belongs to two vertices, while 
each external one to one, so that the sum over all vertices (w) of the number of legs for each vertex L„ is equal to 
L = ^^ Lv = 2iV + E. Adding the contribution from the integrals in the internal lines, this diagram scales like 
X-(i+a) L/2+2N ^ yi^i-a) N~{i+a)E/2^ rpj^jg ^Yiows that cach Order in t^ contributes a factor A^"" ~ Sf'^ (Ref. |2|). 
To get the same diagram in momentum space one has to integrate over E — 1 external xi\ and t, getting a factor 
^2_E-2^ For example, a momentum-space diagram of order t±-'^ for the inverse self-energy scales like £\\^ +i)(i-q)^ 
This is correct provided no short-distance dive rgen ces occur in the integration of diagrams, i. e. if the integrals do 
not depend on the short-distance cutoff a of Eq. (|A9|). A short-distance divergence would introduce a negative power 
of a, which has to be compensated by a positive power of A in order to have the correct dimensions (powers of a lenstii 
scale). This is what happens, e. g. in diagrams 7 and S in Fig. nlc, for a > a2p, i. e. in the two-particle regimell3. 
In diagram 7, if one assigns to the external lines (1, 2) the index r = -t-1, and to the internal lines 3, 4, 5 t he in dices 
r = -1-1,-1,-1, respectively, and inserts the expressions for the two two-particle vertices taken from Eq. ( [A14| ), one 
obtains for r(xi — X2) a contribution of the form (we don't consider the dependence on the _L coordinate here) 



a'"i±' /n^'x, (|xi-X3||x4-X5| 



-1-Q 



4=3 



|Xi -X5I IX3-X4I 



-B 



xe 



xi -X4I 1x3 -X5I 

i(arg(xi-X3)-arg(x4-X5)) 

|X2 -X5I |X3 -X4I 
|X2 -X4I |X3 -X5I 
i(arg(x3-X2)-arg(x5-X4)) 



(A16) 



- 1 
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According to the scaling analysis carried out above, the contribution Eq. ( A16| ) should behave like a^" tj_^ |xi— X2 



|2-4a 



(notice that this expression correctly has the dimensions of an inverse length) . This behavior is correct by assuming 
that the integral does not depend on a in the a -^ limit. However, this is not the case for B > I, for which the integral 
diverges at small distances, as one can readily verify. Thus, for B > 1 the integral gives an a-dependent contribution 
^2-2B ^ which must be balanced by an additional contribution ex |xi — X2 p'^~^ in order to have the correct dimensions. 
Thus, for i? > 1, corresponding to a > Q;2p, the contribution Eq. ( [Al(j ) goes like t±_^ q2-2B+4q |^^ _ ^^\2B-ia ^ -^ ^ 
a stronger divergence. This produces the two-particle exponent obtained in Ref. M. 

APPENDIX B: RESULTS OF THE LOWEST-ORDER APPROXIMATION 

In this section, we summarize some results of the LO approximation Eq. (O) introduced by WenEll. Within this 
approximation, the introduction of t_L modifies the denominator of the Green's function by a term t^c_\_. The Green's 



function for the LM Eq. (C8) can be readily analytically continued to the complex plane (we set the constant ga to 



1 for simplicity, and take r — +1): 

g {k\\,z = iuj) = : ^- = ; . (Bl) 

lUJ — fc|| Z — K|| 

This expression is analytic for z on the real axis and — fcy < z < A;|| , which is the reason why the LM spectral function 
is zero in this region. The denominator of Eq. (0) becomes 

g"(fc|| ,z)-'- c^t^ = (T^T^ - '^'^ ■ (B2) 



The zero of Eq. (B2) gives a true pole whenever it occurs within the region of analicity. For example, the FS is given 



by the points k\\ , cj^ where 

^°(fc||,z = 0)-i=ciii, (B3) 



1. e. 



-fc||/|fc|||" = c_Lt_L^fcj| =-signcj^ |c±t±r/^^~"' . (B4) 

By including a finite value for the energy z, one can easily see that, whenever c± 7^ 0, Eq. (|B2|) is analytic in a 
neighborhood of this point, i. e., the solution is a true pole (cf. Ref. Bq). By differentiating Eq. (p32|) with respect to 



z and replacing the solution Eq. (B4), one obtains the inverse of the residuum, i. e. of the weight Z, for this pole. 
The resuh is Z = (c_Lt_L)"/^^""\ and is plotted in Fig. |. 



Close to the FS, one can look for a zero of Eq. (B2) of the form z = x k\\. This gives 

/c|| =y(x) fc||^(ci) (B5) 

with 

The solution is real, and thus it gives a pole, for each — 1 < x < 1. In this region, y{x) takes all the values 
< y{x) < 00, i. e., for each y > there is always a solution x. This means that for any point (fcy, cj_) in the Brillouin 
zone with kn c± < one always has a pole at a given frequency. The weight Z of the pole is readily evaluated as 

Z = Ic^t^r/ti-") ((1 - a;)-"/2(i + ^)i-"/2) '^^'""^ (1 + (1 - a)x)-^ , (B7) 

which vanishes only at the border of the region, x —f ±1. Obviously , the above discussion only holds for a < 1. 



The fact that there is always a true pole for any A:||Cj^ < can be also seen directly from Eq. (B2). For given c±tj_ 
(say > 0), and q = — fc|| > 0, the function 

(B8) 



(g2 - z2)«/2 

vanishes for z = —q and diverges for z —* q^ . Between —q and +q, it i s an increasing function of z. Thus, for any 
t±c±, there is always a value of z within the analytic region of Eq. (p3), giving a zero. In practice, for small c± the 
pole starts to build close to the left nonanalicity, while for increasing c_i it approaches the right singularity. This pole 
can be seen in Fig. ||. 
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APPENDIX C: EVALUATION OF THE DRESSED HOPPING 

In this section, we evaluate the long-distance behavior of the-, dressed hopping in real space, which we need in 
Eq. (^. Its diagrammatic equation is given in Fig. gb and readscHl 

n{k, Ci) = t^cj_ + t^c^ g(k, Ci) t^c^ = t^c^ (1 - t^c^ r(k))-i , (CI) 

where we have used the Dyson equation Eq. (0). At the lowest order, F scales as F(k) ^ £||"^^, and thus, from 
Eq. (|Cl|), 7j_(k, cj^) formally goes like i^n"'^ " for small energies and fixed tj^cj^. As discussed in the Introduction (cf. 
also Sec. |A|), every order in t± in the perturbation expansion carries along a term which scales like f||"~ and thus 
higher orders in ij_ are stronger and stronger divergent. However, due to the scaling of Tl , replacing the bare hopping 
ti_ in the perturbation expansion with the dressed T±_ cancels this power-law divergence. Thus, the correct starting 
point to study the low-energy region is to carry out a skeleton expansion in T± and remove all self-energy insertions. 
In our case, this corresponds to replacing the diagrammatic series of Fig. |l|d with the one of Fig. ^ja. Although the 
power-law singularities have disappeared in this way, logarithmic divergences are still present in this expansion as 



discussed in Sec. [II. These divergences can, however, be jesummed, in the same spirit as it was done for the parquet 
series by Dzyaloshinskii and by Nozieres and coworkersuilcS, as discussed in Sec. Q 

The behavior of 71 discussed above holds for nonzero cj_. In Eq. (||) we need the a;_L = hopping, i. e. we have 
to integrate over cj^, including cj^ = 0. One should, thus, treat this integration point with due care. We first Fourier 
transform in the _L direction, obtaining 

ri(k,xi - 0) = I dci_ i?(ci)ri(k,ci) , (C2) 

where T>{c±) is the density of states-fpt the out-of chain energy. In the Z? — > oo limit and for a cubic lattice with 
nearest-neighbor hopping this readsES^La 

V{c^) = -^e-^^'l^ . (C3) 



The integral Eq. (02) can be readily evaluated for small energies, where the quantity e = [tj^r(k)] is small. By 



inserting Eq. (CI), collecting e and summing and subtracting i_L yields 



r^(k,x^ = 0) = -e j dc^V{c^)(t^ + ^-^^ =-£t^[l + 0{e loge)] ^ -^ , (C4) 



where the logs contribution is given by the cj^ = point. It is clear that the asymptotic result Eq. (C4) does not 
depend on the specific form of the density of states I?(cj^), as long as it is regular at c_l = and normalized. 

We now carry out the Fourier transform in the || direction. For the sake of definiteness, we consider here the right- 
moving (r = +1) component. All results for r = — 1 are simply obtained by changing the sign of the || coordinate, i. 
e. x\\ or fell. As a first go , we evaluate Tl'^, the LO approximation for 7 j_ , i . e., we use F = Q'^ . The full dressed 71 



will be evaluated in Sec. C 1. The LL Green's function is given in Eq. (A12). Its Fourier transform is given by 



g"{k)= fd^xe-'''-^g°{x\0), (C5) 

where we have identified^ k = (fc|| , — cj). We now introduce the angles of the two vectors x and k with the "x" axis, 
i. e. (f> — argx • v, and 6 — argk • v, and v = (1, i) as in Sec. ^. Eq. ( |C5| ) becomes 



'ta 



^O(k) ^ ^^ / ^2^ g-»|x||k|co.(e-« g-»0 |^|-l-a (C6) 

27r 



Going over to circular coordinates and transforming (p' = (p — 9, and s = |k||x| yields 

poo p27T 



g\k) = |k|"-i e"'^ / s~"ds / d0' e-'('^'+^™^'^') , (C7) 



StJ" J|k|a Jo 

where we don't care about the specific form of the cutoff at s < |k|a, as it can be taken to zero in the low-energy 
limit (|k| is always limited by teff <C 1/a)- The last integral over 0' gives —27riJi{s), with Ji{s) a Bessel function. 
Integrating over s, Eq. ( |C7| ) gives 
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where 



k • V k\\ — iLO 



_ T{l^a/2) 
5a-2ar(l + a/2) ■ ^ ^ 



Like in Eq. (C8), we will from now on indicate with "~" expressions valid in the asymptotic limit. However, whenever 
this will become clear, we will switch back to "=" . 

To evaluate Tl°(x, 0), we first insert Eq. ( |C8| ) in Eq. ( |C4| ) (remember, here we use F — t/O), and then transform 



back into x coordinates. Thus, 

Ti°(x, 0) « - y ^ e*'^-'' g\\^)-^ (CIO) 



with the same conventions as above, and with g = |k|. Transforming q\x.\ = s and integrating over 0, yields 

I |Q-3 10 j-oo 

Ti"(x,0)«^ — / s'-^ ds 2^1 Ji{s) . (Cll) 

47r-^ g^ a" Jo 

In principle, the last integral does not converge at large s. However, it can be regularized by inserting a convergence 
factor e"'" = e"''!^!!''! with ^ - ^ — , physically due by the fact that the behavior TL(k, c_l) ~ -r(k)-i [ Eq. (|c|)] 

is cutoffed at |k| ^ teff- The convergence factor fi can then be safely taken to zero, since the result of the integral 
does not depend on /i for small /i. In this way, one obtains 

r^°(x,0)«^^^^|xr3e»^ (C12) 

with (p — argx • v. 

1. Fully dressed function 

We now carry out the same Fourier transforms with the renormalized function F, i. e. with ( Eq. (|6|) with m = {]) 

F(x) = g,(x|0) = gO(x|0) .Fo[alog(|x|ie//)] , (C13) 



with the renormalization function !Fq{1) given as a power expansion (cf. Eq. (Fl)) 



•^o(0-E/"^"- ('^l^) 



n=0 



The Fourier transform of F can be carried out as in Eq. (C7), and the integral over 0' gives the same result, as Tq 
only depends on the modulus of x. Thus, we are left with 

r(k) = ^— |kr-i e-*M s-"ds[-2^zJi(s)]V/„a" logs + log^ . (C15) 

Since !Fo is, in general, a complicated function, and its coefficients /„ very general, the Fourier transform can be only 
carried up to the leading logarithmic behavior, which, as discussed in Sec. Q, amounts to considering I = alog(|x|te//) 



of order 1 but a small. In this way, we can neglect the logs within braces in Eq. (C15) and the effect of the 
renormalization function Tq becomes merely multiplicative, provided one replaces |x| with |k|~"^ in its argument. We 
thus obtain 

r(k) = -^^o(alog^)=e°(k)^„(^alog^) ' ^^^^^ 
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where we have replaced the coefficient g^ with its a ^- Hmit (7a=o = l, consistently with the leading- lo garit hmic 
approach. One can also verify a posteriori that inserting the asymptotic result for J-q{1) (^ e ) in Eq. (CIS) one 
indeed obtains Eq. ( |C16| ) for small a. 

We now need 71(x, 0) in real space, i. e. the Fourier transform of — r(k)^^. To express — r(k)^^ we need the 
reciprocal function of J-q in terms of its power-series coefficients /„: 






^o(0 



(C17) 



n=0 



where /„ can be determined from all the /„ with m < n. Again, this function does not depend on angles, and we 
can proceed as for Eq. (Cll), yielding 



71(x,0) 



IQ— 3 i4> /"oo 



47r2 ga a" Jq 



'" ds 27n Ms) J2 /" «" [logd^lW/) - logs]' 



(C18) 



n=0 



The procedure is now slightly more complicated than for Eq. (C16), since we have to con sider terms at the ffist order 
in logs. The reason is that, if we neglect completely the logs term, the integral in Eq. ( Cl^ ) is of order a: 



2~a 



Ji{s) ds = ga- 



2a 



(C19) 



On the other hand, expanding the [• • •]" power on the r.h.s. of Eq. ( |C18| ), and keeping the first term in logs yields 
a result of order 1: 



f°° _ d 

/ s^""logs Ji(s) ds = -— ga_2 « -2 . 



(C20) 



The first integral, thus, gives a contribution 2 /„ q;"+^ [log(|x|ie//)]" to the n-th term of the series in Eq. (C18) 



while the second gives 2 n /„ a" [logdxlie//)]" . Both terms are of the same order within the leading-logarithmic 
approach and must be taken into account. We thus obtain 



Ti(x,0) 



n a" 



|x|"-3 g»0 



J^o(alog(|x|ie//)) +^0 (alog(|x|ie//)) 



(C21) 



since XI^^o /" (^" + ^ '""^) = -^0(0 + ^o (0 (here, To (l) == ^^o(0)- Eq. ( |C2l| ) is the final resuh of this Section, 
which we need to insert in Eq. (0). 



2. First 1/D corrections: irrelevance of _L— nonlocal dressed hopping 

In the Z? — > oo limit, only the local effective hopping T±{x.,x± — 0) is needed in the diagrams of Fig. 0, as _L- 
nonlocal contribution vanish in this limitt^. In order to study the contribution of finite-Z? corrections, we consider 
the _L-nonlocal contributions to T± , given by 



T^{k,x^^O) ^ I dc^ V^^ic^) rL(k,Ci) , 



(C22) 



where, since T± depends on k± only through cj_, we have introduced the "generalized density of states" (here, we use 
R instead of x± ) 



Vr{c^) 



D' 



n 4^ e^'^^^ \ 5\c^ 



D' 



2tt 



/D 



= ^ cos kd 



rf=i 



Following Refs. 39, we now introduce the Fourier representation of the S function, obtaining 



(C23) 



(C24) 
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where the integral /(s, R) is given by 
Hs,R)^ [H ^ e'k,R,^^^scosk,/^ ^Y[ [ ^ e''^''^ (l-^cosk-^ cosH + 0{D'-^^')] , (C25) 



and we have expanded in powers of l/^/ly. The last integral gives at the leading order 

1 



D' 


for Rd = 


= 


7W 


for Rd - 


= 1 


D' 


for Rd = 


= 2 



(C26) 



and, in general, a term of order (s/v-D')" f^^' -^d — "• Inserting these results in Eq. (C24), one obtains at the leading 
order in Xjyiy 






ID' 



Rd 



where the an are coeffici ents obtained from Eq. (|C25|). For example, from Eq. ( ]C2q ) oq — 1, ai = —i, 02 
The powers of s in Eq. (C27) can be replaced with derivatives with respect to cj_, yielding 



^i?(c_L) = !!"«<* 



fjy dci_ 



2?(c±) , 



(C27) 
-1/2. 

(C28) 



where the usual density of states is given in Eq. (|C3|). We can now insert Eq. (CI) and Eq. ( |C28[) , in Eq. (C22). 
Since J dc± I , ) I'(cj^) = for n > 1, for the nonlocal T^ we can subtract a cj^-independent term from 7j_(k, cj^), 
and write 



71(k, xj^ 7^ 0) ex / dc±^ 



ri(k,ci) 



r(k) 



1 



/Jj dcj_ 



T^{c±) , 



(C29) 



where x± = ^^ x±d- The term within brackets in Eq. ( |C29[ ) goes like r(k)^^ for large r(k), and the same holds for 
the integral (the fact that the coefhcient of r(k)~^ diverges at c_l = might, at most, give a log correction). Thus, 
71(k, a;j_ y^ 0) vanishes at least like r(k)~^ for low energies, i. e. faster than T±{'k,x± — 0). Diagrams containing 
_L-nonlocal T± contributions are thus irrelevant in the renormalization-group sense. 



APPENDIX D: INTEGRATION OVER CENTER OF MASS COORDINATES 

In this section, we prove Eq. (R), for the integration over the center of mass coordinate ym+i- In order to 
simplify the notation, we introduce the shorthand C(0, • • • ,71) = ^^(yoi ' ' ' lYnlyoj ' ' ' ^yn), for the cumulants, and 
G(0, •••,n) = CyO(y[,, • • • ,yJJyo, • • • ,y„)j-ipr the disconnected LL Green's function. Notice that in these Green's 
functions the implicit a and r variablesEHI are pairwise equal. More precisely, ak and r^ , associated with y^, are 
equal to crj. and rj,, associated with yj.. The reason is that a Tl (or t±) line does not change neither a nor r, see 

Fig. ||. In addition, we define x^ = yj, - y^, F{n) = 27r|x„p Gin), Ij^n = "log j^ {^rj,r„ + ^Sr,-r„), and for 

the integration J (i^y„ we use the notation J . 

The proof proceeds in two steps. We first show that the term on the r.h.s. of Eq. (|^) is also given by an integral 
of disconnected Green's function, namely 



[G{0,---,n)-G{0,---,n-l) G{n)] 
G(0, •••,n 



(Dl) 



n-l 

1) F{n) J2 kn , 

J=0 



where, for convenience, we have renamed rn — > n— 1. Then, in Sec. p 3| , we prove the step from Eq. (Dl) to Eq. (M) 
by induction. 
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Disconnected Green's function 



To show the first part, we write the n+ 1-particle correlation function Eq. ( All ) by using Eq. ( |A2| ) in the following 
form 



n-l 



G(0, • • • , n) = 0(0, • • • , n - 1) G(n) x J| 



^n,r„(y» - yn)fi.,r„(y^ ~ y«) 
Pr,.r„iyi - yn)-Pn,r„(y- - y!i) 



where we have used the fact that rk = r'^- We thus have 

[G(0, • • • , 7i) - G(0, • • • , n - 1) G{n)] = G(0, • • • , n - 1) G(n) / I{n ~ 1) 
with the argument of the integral 



/(n-1) 



n 

4=0 



-Pn,r„(y» - y'n)Pri,rAy'r ^ Yn) 

^n,r„(y4 - yn)i^n,r„(y- - yJi) 



1 



(D2) 



(D3) 



(D4) 



The integral in Eq. (|D3|) is restricted to the region -IJ- n, where x„ is smaller than all other distances, which are 
the arguments of the P's, in Eq. (D4). For this reason, we can expand I{n — 1) in powers of x„ = y^ — yn- The 
zeroth order of this expansion is zero, as I{n — 1) = for x„ = 0. The first order I{n — 1)^^^ gives 



/(n-l)W = 5: 



i=0 



V-Pn,r„(y» -y«) 

Pri,rAyi -yn) 



V-Pn,rJy--yn) 
^n,r„(y- -y«) 



(D5) 



where the V is considered as applied to the argument of the function. Having in mind to integrate this expression 
over y„, one would be tempted to carry out a shift in coordinates yn — > yn + x^ in the second term within brackets 
in Eq. ( |D5|) , thus obtaining zero. This shift, however, has to be carried out with some care, since the logarithmic 

gradients p '''■''y^^ ^y go like l/|yn| for large |y„|. Therefore, the integral of each separate term in Eq. (Dq) 

does not converge, i. e. the shift is not allowed without further prescriptions. However, this hol ds i f one uses the 



zero-temperature form Eq. ( |A9| ). On the other hand, the finite-temperature prescription Eq. (A7) introduces a 
cutoff for values of each of the arguments in the VP/P in Eq. (D5) of the order of 1/T, making the separate integrals 
absolutely convergent and allowing for the coordinate shift. 

We thus need to expand /(n — 1) up to the second order in x„: 



1 + /(n - 1) « n (1 - < ^/^(y^) + < < ^M.-(yO) n (i - < P^^'^y'j) + < < P^'y^S'i) 

3=0 



i=0 



(D6) 



where a sum over repeated indices fi, //', v, v' is understood, and where we have introduced the notations y^ = y-i — yn, 

and y- = y- — yn- Moreover, x^ is the /i component of the vector x„, Pfi{y) = p, \ , and P^,u{y) = 

' ^ 2 P(^) — Moreover, we have omitted the r^ indices in the functions P, since they are fixed by their y arguments 

[i. e., P{yi - y„) = Pn^r^iyi - yn)]- 



Expanding the denominator of Eq. (D6) we obtain 



1 + /(n - 1) « n (1 - < ^A'(y') + < < P^^Af^)) 

i,j=0 



(D7) 



x(^i + x^ p^,(y;)-x^ K P^',.'(y;) + x^ x^p^-(y;.)P.'(y;-) 

Collecting powers of x^j, we obtain the second-order term /(?i — 1)*^^^: 

i{n - 1)(2) = x^ xj; {Y^ [p^Afi) - P.Af'i) + PM) PM)] 

i 

+ J2[p,{y.) PAyj) + PM) PAy'j)] -E^M(y') ^--(y^)} 



i>j 
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i 



(D8) 



With the integration over y„ in mind, and with the same arguments abo ut convergence as for Eq. (D5), we can carry 
out a coordinate shift of y„ in some of the terms of the sum Eq. (D8). First of all, we shift y„ —t y„ + x^ in the 
Pti,v{Y'i) term, so that it becomes P^,,y(yi) and it cancels the first P^,u term. Next, we transform the first-derivative 
term in the first sum in Eq. (D8) in the following way 



{PM) - P^^{yi)) PAfd 

= I PM) PAfd + 1 PM) P^y'i) - \ P^{y^) pM) - I ^^(y.) ^.(yD 

^ I PM) PMi) + \ PM) PM) - \ PM^) PMi) - \ PM'd PM) 

= \ [PM'i) - PM'^)] [PM)-Pu{y'i)] , 



(D9) 
(DIO) 



where Eq. (DIO) is obtained by shifting yn ^ yn + x.^ in the second term and by exchanging /i and v in the fourth 
term of Eq. (D9) (which is allowed, as Eq. (D^) is symmetric in /i, t^). Inserting the result in Eq. (D8), and 
factorizing in the same way the terms in the last sum, we finally get 



i{n - 1)^^' -^ \ < K E [^^'(y^) - PM'i)] [PAy,) - PAvj)] 
= l\j2< lPM^)-PM'i)]\ , 



(Dll) 



where "^" means that it is equal but for a shift of the integration variable y„ in some of the summands. 

We now need the logarithmic gradients P^{ya ^ yt) - If the points a and b correspond to two electrons on opposite 
side of the FS, i. e., Va = -n, then from Eqs. ([A^, A4pi| ), P°(y) = lyl""^, and 



PM) - %^ 



-B 



-B- 



(D12) 



where we have set y = y^ — yf,, and introduced a superscript symbol or 1 to P^, depending on whether Tq — — r^ or 
^a = fb, respectively. In the second case, Ta = rf, = r, we can write P^{y) — c |y|~^~" y • v, where c is a constant, and 
the two-component vector v = (1, — i r), slightly different from the one defined in Sec. 1^. Differentiating, we obtain 



PM) - -(2 + «)fii + ^ = n2 [(y • v*K - (2 + «)yn , 
|y| y-v |yr 



(D13) 



as |yp = (yv) (y-v*). 

There are thus three types of integrals to be carried out in Eq. (D3) with the second-order term Eq. ( Dll ). First, 
for the case that r„ = — r^ = — ''j, we need an integral of the form 



pM^ - y") ^"(y^ - y«) = ^' / ^^y 



■xP|y|2 



= B^TT 5i,y log 



R 



max(|y, -yj|,|x„ 



(D14) 



where in the intermediate step we have transformed yn = y + y^, and y^ — y^ = x, and used the result Eq. (D21). 
Here, we have introduced a large-distance cutoff R ex 1/T, from which, eventually, the final result does not depend. 
The max in the logarithm practically only applies when j = j, as jx^j is always the smallest distance. In this case, 
the result is obtained by keeping in mind that |x„| is the short-distance cutoff, and by applying Eq. (D20). For 
fn = Ti ^ Tj ~ r, we need 



I pM^-y^)PM,-yn) 

J 11 



(D15) 
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1 



|yP|y-xP 



{[(y - x) • V*] v^ - (2 + a) (y" - x'^)} [(y • v*) v'' - (2 + a) y-^] = 



v*M v'^ - (2 + a) 6^^ 
= TT (2 + a) a Sf_,^ log 



V- v'' - (2 + a) S,,, I d^y^ ^ 



y-xp |yp 



R 



max(|yi -yj|,|x,, 



again using Eq. (D21) and the fact that v*'^ v*'^ — 0, and v^ v**^ + v*^ v" = 26^,^. Finally, for r„ 



have 



P"^{y^-yn)Pl{yj-yn) = {-B) 

Tr[B{2 + a) (5^, - B v*^ v^ log 



(2 + a) 5^'^ 



TT (^u,!/' log 



R 



R 



max(|y,, -yj|,|x„ 



TT (1 + a) B 5^,, loi 



max(|y,; -yj|,|x„| 
R 



-r,, we 



(D16) 



max(|yi-yj|,|x„ 



where, in the last step we have symmetrized with respe ct to /i and v. 
We can thus use these results to integrate Eq. ( Dll ) and obtain 



r ^^^ Iv' 

//(n-l)(^)«5<5,.xj:x^{^ log^ 



yil|y»-y 



i^j=0 



vjm 



fj [((2 + a) a(5^.,^„+B2 Sr„-r^ 



I I 
+ ^2 ((2 + a)a<5,.,,„+B2^,.,_,J log^j, 



(D17) 



i=0 



where we have considered the case i = j separately, and used the fact that |x„| is (niuchEj) smaller than all other 
distances in the region ij. n, and, thus, it can be neglected whenever it appears summed to other distances as the 
argument of a logarithm. Notice t hat th e large-distance cutoff R cancels out, as anticipated. 

Consider now the terms in Eq. (D17) with i =/= j. These give logarithmic contributions of the form 



yj+Xi -y^ yj -yj -Xj 

log! n — — 

y* -yj yi + xj -y^ -x^ 



(D18) 



For the sake of definiteness, let's take i > j in Eq. (D18), so that, in the relevant region ])- n, Xi is smaller than 
all other differences i n th e arguments of the logarithm and thus can be set to zero. In this way, numerator and 
denominator in Eq. (D18) cancel and the result is zero. This means that the terms wit h i ^ j in Eq. (D17) do not 
contribute to the leading logarithmic divergence. Thus, the only contribution to Eq. ( D17 ) stems from the second 
term within brackets, which gives 



//(n-l)(2)«2^|x„p^alogf 



(D19) 



where we have taken B^ w 2a/ S, and (2 + a) a ~ 2a, consistently with the leading-log approximation. Inserting 
Eq. ( pl9| ) into Eq. (D2) yields the desire d result Eq. (Dl). • • • well, there is certainly a much faster and elegant way 
to get the rather simple result Eq. (D19)!. 



2. Some logarithmic integrals 



Here, we evaluate the integrals used in Eq. ( D14 ), and following. The first integral is straightforward 

,2 y^y'' c , R 

^ y I^aT = ^ ^°'P log X ' 

'A<|y|<fi |y| ^ 



/ 

JA 



(D20) 



where i? is a large-distance and A a short-distance cutoff for |y|, which are needed due to the logarithmic divergences 
of the integral. We next prove that 
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d y T 



|y|<« 



i2 j:a2 - "" ^'-13 logo 



(D21) 



where « means at the leadmg order in log(i?/|x|). The integral converges at short distances, so there is no need for 
a short-distance cutoff, and diverges logarithmically at large distances. We can split the integral in two regions: (i) 
|x|A^ < |y| < R, and (ii) |y| < |x|A^, with N large but much smaller than i?/|x|, so that logiV can be neglected. In 

region (i), the integrand can be safely approximated by ^ u , whose integral, taken from Eq. ( D2Cl| ), gives tt Saf3 log -A. 
In region (ii), the only length scale left is |x|, since the integral converges at short distances, and thus there is no 
logarithmic contribution from this region, and Eq. (D21) is proven. 



3. Integration of cumulants 



We have thus proven Eq. (Dl), an equation similar to Eq. (|7|), but with disconnected Green's functions instead 
of cumulants. We now prove by induction the same thing with cumulants. Induction is the best way to do it, as 
cumulants themselves can be written by induction in terms of disconnected Green's functions. A n-particle cumulant 
consist of the sum of the n-particle disconnected Green's function plus an appropriate sum of products of /c-particle 
Green's functions with k < n (Ref. 36). However, we can show that in our problem, we can restrict to the so called 
"paired" contributions to the cumulants, i. e. we can throw away all those terms in the sum in which, for any fc, 
the coordinates y^ and yfc + x^ do not belong to the same Green's function. The fact that these terms ("unpaired 
terms", see Fig. nO) can be neglected is shown in Sec. D4. 

Let us write in the shorthand form 



C(0,---,n)«C(0, 



1) F{n) 



n-l 
3=0 



U 



(D22) 



which coincides with Eq. (0) for n = m + 1. The induction procedure consists in proving (i) that Eq. ( D22 ) holds 
for n = 1, and (ii) that in the hypothesis that Eq. (D22) holds for all n < tti, it also holds for n = to + 1 . 

For n = 1, C(0, 1) is equal to G'(0, 1) — G(0)G(1) plus u npai red terms. Since, as discussed abov e, u npaired terms 
can be neglected, for n ^ 1 E q. (D22|) coincides with Eq. (Dl), which we have just shown in Sec. D 1. 



We now assume Eq. (D22) to be valid for all n < m. Let us first introduce the definition of a cumulant in terms 
of connected Green's functions. 



C(0, •• -,71)- G(0,- ••,«)- 



E 



C{Pi)---C{Pn,) 



(D23) 



where we have already left out unpaired terms. In Eq. (D23), the P^ are subsets of the set of integers {0, • • • ,n}. 

y ) goes over all inequivalent partitions 



{Pi, • • • , Pnp} is a partition with Np terms of this set, and the sum 

p(o, •••,") 
with Np > 2 of this set. Equivalent partitions are the ones which can be set equal by a permutation. Introducing 
Eq. (D23) in the result for the disconnected Green's functions Eq. (Dl) (with n ^ rn + 1), one obtains 



rn+l 



C'(0,---,to + 1)-G(0,---,to)G(to + 1)+ J2 C{Pi)---C{Pnp) 

P(0,--,m+l) 



(D24) 



G{Q,---,m)F{m+l)Y,h.m+i 



In Eq. (D24), the sum over the partitions of the set of integers 0, • • • , to + 1 can be further splitted in the following 
way 

Np 

Y, C{P,)---C{Pnp)= Y. Y.^{Pi)---C{Pk,m+l)---C{PNp) 

P(0,---,m+l) P(0,---,m) k=l 

+ E C{Pi)---C{PNp)C{m+l) + C{fd,---,m)C{m + l) , (D25) 

P(0,---,m) 
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i. e., into the sum over the partitions of the integers 0, • • • , m with t he el ement m + 1 either appended in all subsets 
of th e par tition, or taken alone. Upon applying the definition Eq. ( D23 ) with n = to to the last term on the r.h.s., 
Eq. (D25) becomes 



Np 



Y, C{Pi)---C{Pnp)^ Y. I]C'(Pi)---C(Pfe,TO+l)...C(Pw,) + G(0,---,TO)C(TO + l), (D26) 

P(0,---,m+l) P(0,---.rn) fc=l 



which, inserted into the l.h.s. of Eq. ( D24 ), cancels the second term within brackets, giving 



C(0,---,TO + 



m+l 



. Np 

1)^-/ E Y.^{Pi)---C{Pk,m+l)---C{PNp) 



(D27) 



P(0,---,m) fe=l 



+G(0, • • • , to) F{m + l)Y ^J^™+i 



The integral J . in the first term on the r.h.s. of Eq. (D27) can be evaluated by using the induction hypothesis 
Eq. (D22) with n < to, as C{Pk,m + 1) is a cumulant with less than m + l particles. We thus obtain for this term 



Np 

E E^(^0' 

P{0,---,m,) fc=l 



C{Pk)F{m +l)Ykm+i 



jePk 



■C{P, 



Np 



(D28) 



■ Y. C{Pi)---C{Pk)---C{PNp)F{m + l)Ykm+i 

P{0, — ,m) j=0 



^Np 



since X]fe=i Siep — T^T=o- Inserting the last result in Eq. ( D27 ), and using again the definition Eq. ( D23 ) yields 



the desired result, i. e., Eq. (D22) with n = to. + 1 



4. Irrelevance of "unpaired" terms 

We want to show that the "unpaired terms" in a cumulant do not contribute to the leading logarithmic divergences 
in any of the terms of the sum Eq. (H). A ?7i + 1-particlc cumulant is the sum of products of n-particles Green's 
functions with n < to, + 1. For "unpaired terms" we mean those terms in the sum for which some paired variables (i. 
e. Yk and y'^. = y^ + x^) do not belong to the same Green's function. For example, for to, = 1 



G2(0, 1|0', 1') = G°(0, 1|0', 1') - G°(0|0') G"(l|l') + G°(0|l') G°{l\Q') 



(D29) 



the last term on the r.h.s. is unpaired, while the first two are paired. The first two terms, when inserted in Eq. (^ give 
a log^(|xo|te//) contribution, as discussed in Sec. III. The contribution to F of the last term can be best understood 
diagrammatically (see Fig. |lO| ). Its contribution, written in momentum space, is proportional to 



I (f\i G"(k) rL(k,2;_L = 0) G"(k) e*-^° ^ - f (fk 



GO(k) 



^o(alog^) 



zk-XQ 



GO(xo) 



J^o[alog(ie//|xo| 



(D30) 



i. e., it does not give additional logarithmic terms to F, contrary to the contribution from the paired terms. The 
same thing happens at higher order, namely, while the paired terms of a to + 1-particle cumulant inserted in Eq. (jj) 
give a correction of order log '"(|xo|te//) to F, the unpaired terms give smaller powers of the logarithm and can thus 
be neglected at the leading logarithmic order. 

APPENDIX E: RELEVANT INTEGRATION REGION 

In this section, we show that the leading contribution in log(|xo|ie//) to each of the integrals in the series Eq. (|j), 
restricted to |xm| < |x„j_i| < ••• < |xi|, can be further restricted (A) to the subregion |xi| < |xo|, and (B) to 
|xpl < |yg + ^i^q ^ Yq' ^ e2Xg' I for each p> q,q' , q ^ q' , and ei = 0, 1. This relevant region, which we call "0 -|| to" , 
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0' 




r'^ /I = ,, I ;i + i'>^- - - _j 



0' 0' 0' 



FIG. 10. Splitting of a cumulant contribution into "paired" and "unpaired" terms. A cumulant is indicated by the black 
dot and is obtained as a sum of disconnected Green's fu nctio ns, represented by black squares. The last diagram on the r.h.s. 



is an "unpaired" one, according to the definition of Sec. D4, and does not contribute to the leading logarithmic divergences, 
as shown in that section. 

is the one where the modulus of the relative coordinate Xp with a given index p (p = 0, • • • , m) is smaller than the 
distance between any two different points with indices q, q' smaller or equal than p. The remaining regions do not 
contribute to the leading logarithmic divergence of the integral. This is a crucial point in giving the simple expression 
Eq. (I). 



Let us start with m — 1. We have already seen in Sec. Ill that the contribution to Eq. M) from the region JJ- 1 
gives a log term, and, for general m one has a log ™ contribution. Consider now the integration region |xo| < |xi| 
in the term m = 1 in Eq. (0), which violates (A). The integration over the "center of mass" coordinate yi can be 
replaced with an integration over yo, since the integrand depends on the difference between the two. As a consequence. 



one can just take over the result Eq. (Dl) with n = 1, and interchange the labels 1 and 0. This is correct because 
now |xo| is smaller than |xi|. One thus obtains 

/ [G(0, 1) - G(0) G(l)] ex 2 a TT G(0) G(l) |xop log p\ . (El) 



If one now inserts the expression Eq. (012) for the LO 71, and integrates xi from |xi| = to |xi| = |xo|, the result 



is proportional to a^ G(0) |xoP /, i^, , ^^ log |^ — a^G(O) x 0(1), where 0(1) is a term of order unity, i. e. 
without logarithmic contribution. For a generic term of the log expansion, Eq. ( |C2l| ), of the completely dressed T± 
one has a similar result, namely after integration over xi one has no additional log contribution, while one gets a 
term a^, i. e. one "looses" two logarithms from integrating in that region. 

Taking now m > 1, one first integrates the variables y2, X2, • • • , ym,^m by using Eq. (H). However, this integration 
simply renormalizes the two-particle cumulant Qc by a factor of order 1 in the leading log, i. e. by a sum of powers 
of (alog |x|te//). Now one can proceed integrating over yi,Xi. By the same argument as above, it is straightforward 
to show that integration from the region |xo| < |xi| does not contribute additional logarithms, while it gives a term 
a^ , and can thus be neglected. We have thus proven (A), i. e. that the leading divergent contribution to each term 
of the expansion Eq. (0) comes from the region |xo| > |xi|. 

To show the second part (B) of the statement, we should firs t understand ho w a l ogarithmic contribution to Eq. (H) 



comes out. Let us consider the integration over y™ of Sec. D 1. I{m~ 1) ( Eq. ( P4D ), which is the only y„j-dependent 
part of G(0, • • • , m), can be written in the generic form 

/(m-l)=[]7'.(y™-r,), (E2) 



where the Vi are functions with an integrable singularity in 0, [whenever the exponents 1 + a, and B ( Eq. ( A6)) are 
smaller than 2]. Since the singularities are integrable, there is no divergent (power law or logarithmic) contribution 
from integration of y™ in the neighborhood of the points r^. For the sake of definiteness, let us suppose that ri and 
Y2 are the two nearest points among the r^, and call A = |r2 — ri| their distance. Then, one can consider the circle 
-R<A of radius A^A around ri, with N some number of the order 1 smaller than the relevant logarithmic scale. This 
region contains ri and r2 but none of the other r^ points (In case there are other points r^ inside this region, the 
following argument does not change, provided their distance from ri and r2 is neither much larger nor much smaller 
than A). Inside i?<A, there is just one characteristic length scale A, since there is no need for a short-distance scale 
due to the convergence of the integral. Thus, by simple dimensional analysis, one obtains for the integration in this 
region /„ d^ym li'm — 1) oc A^ with no log contribution since one needs two length scales for a log. A logarithmic 
contribution can only come from integrating y^ in the remaining region i?>A, where more energy scales are available. 



In this region, one can expand in powers of A, when it appears as an argument of the Vi, as we have done in Sec. D 1 
whith A = |x„i|. 

Let us start from the simplest case m — 1. Here, there are four points r^, namely yo, yo ~ xi, y'^ = yo + xo — xi, 
and yo + xo. Since we havella |xi| ^ |xo|, the smallest distance A between two of these points is given by |xi|. 
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As discussed above, the leading logarithmic contribution to Eq. (Q) comes from the region i?>|xi|j where |yi — 
yolJyi — Yoli IYi — YoIj ly'i — YoI > l^il which proves result (B) for m — 1. It is now straightforward to extend this 
argument by induction for any m. Specifically, we first assume that we can restricts to the region where the distances 
\yp - Yql, |yp - y'q\, \y'p - y'q\ (let's call them "|yp - y,| and primed") are larger than |x,„_i|, for p,q<m- 1. Then, 
since |xm_i| > |x„i| (we are restricting to the region 1 | m), it remains to be shown that the region where any one 
of the distances |y„i — yg| and primed is smaller than |x„i| does not contribute to the leading logarithmic divergence. 
Since |xm| is smaller than all other distances |yp — y^j and primed, we can apply the argument above, according to 
which logarithmic contributions from the integral in dPym come from the region outside of circles of radius jx^l from 
any of the points y^ or y' . This proves the statement. 

APPENDIX F: SOLUTION OF THE RECURSIVE EQUATION BY POWER EXPANSION 

In this section, we describe the practical procedure to solve the recursive set of equations Eq. (M) by power 
expansion up to very high order. We also show some results of the corresponding Pade resummation. We expand the 
functions Tm in power of their arguments 



T (II ] — S^ f*-™^ r P 

J m\^: I'm) — / ^ Jij '' 'm 



(Fl) 



We can use two known results, namely (i) JF„i(^,0) — 1, which implies /,^™ = 5ifi, and (ii) fl ■ = Si,o fj, as J-o{l,lo) 
only depends on ^o- 



Inserting Eq. (Fl), and the expansion for the reciprocal function Eq. (C17) in Eq. (hd) yields 



E 47 ^' ^™ = 1 + 2(5+1) / dl' [I + l,n - (to + 1) I'] 



i,j=0 



(F2) 



E f^T'^ (' + imY I" Y.^fp + (P + 1) fp+i] i'^ 



r,s—0 



p^O 



1 + 2(5+1) / " dl' J2 f^s^'^ [/P + (P + 1) /p+l] [(^ + lrnr+' 1''+" - (TO + 1) (? + IraY 1'^+^+'] 



Carrying out the integration and applying the binomial expansion (with the agreement that (J") = for s > r), 
Eq. (P3) becomes 



r+l 



1 + 2(5+1) E /(T+'^ [fp + (P + 1) fp+i] E ^' C^"^'-' 

r,s,p— q—0 



r + l 



q J s + p+1 \qj s +P + 2 



r\ m + 1 



(F3) 



1+E E 

i— J— inax(l,2— i) 






(m+1) 



where we have replaced q = i and r + s+p + 2 — q = j, and introduced the coefficients 

i+j-2 i+j-r-2 



B, 



(m+l) 



2(^ + 1) E E f^r-p-2 [fp + (P + 1) fp+l] 

r=max(i— 1,0) P=0 



r + l\ 1 

i Ji+j — r — 1 



r\ m + 1 
iji + j — r 

(F4) 



Comparison of Eq. (F2) with Eq. (P3) gives the relation between the J*^™^ and the _/'(™+i)^ namely 



„M^ l47^'^ for^+J>2,^>0,J>l 

JiJ S Ij 



0, 



for i — j — 
otherwise ; 



(F5) 



From Eq. (F4) it is not too difficult to prove another restriction on the coefficients /^ , namely /j^''-' = 0. 



p(m) 
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FIG. 11. log[J^o{l)] vs I for S — 2, obtained with a rational Fade interpolation to the expansion of Sec. [F] of the form 
^oil) = Pn{l)/Qkil), with n = fc = 20 (solid line), n = 21, fc = 20 (dashed), and n = 20, fc = 21 (dotted). 
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FIG. 12. d log[Fo{l)]/d I vs I with a logarithmic Fade approximation d log[J^o{l)]/ dl — Pn{l)/Qk(P), with n = fc = 20 (solid 
line), n = 21, fc = 20 (dashed), and n = 20, fc = 21 (dotted, covered by the solid line). 



We have evaluated the coefficients /,- 



/qj up to j = 42, and, consequently, all other coefficients /I™ , up to a 
corresponding high order, by means of an algebraic manipulation program. If one tries to naively sum the series, 
one comes out with apparent divergences already at I of the order of one, which is probably the convergence radius. 
The Pade method is most appropriate for extrapolations beyond the convergence radiua^. Indeed, we find that the 
possible poles are never on the positive real axis, neither close to it, which is the only region where we need J^o{l) to 
be well defined. The simplest Pade interpolation consists in equating the coefficients of the series to the ones coming 
from a rational function Pn{l)/Qk{^), where Pn{l) is a polynomial of order n, and Qk{l) one of order k. In Fig. |ll| 
we have plotted logP„(/)/Qfe(Z), as obtained by this result, with different n,k close to 20. In all three cases, the 
logarithm seems to eventually acquire a constant slope, suggesting an exponential behavior for Tq. However, above 
I ~ 5 the three different interpolations give different results, which signals a failure of the Pade procedure for Z ^ 5. 
A better approximation is achieved by making a rational interpolation to the logarithmic derivative, i. e. to set the 
Ansatz (ilog[jFo(0]/ ^1 — Pnil)/Qk{l)- As one can see from Fig. H^, this logarithmic Pade interpolation now works 
good up to a larger I ^^ 20 and it clearly shows that d \og[J-'o{l)]/ dl ^ 1 for large / (within about 10~* of accuracy), 
i. e. that .Fo(/) ex e'. 

We also want to study the behavior of the interaction vertices, given by the RRC Eq. ^ . To this end, we have 
evaluated their asymptotic behavior, when all internal variables |xfc| are of the same order of magnitude (see the 
discussion in Ref. 46). This is obtained by setting all l^ = I in Eq. (S), or, equivalently, I -^ m I, and Im = I 
in Eq. ([lO|). In Fig. |lj, we have evaluated the logarithmic derivative of Fm{m 1,1). The figure clearly shows that 
Fm{m 1,1) ex 6^™+^^ ', i. e. the associated RRC, Qc of Eq. (ph, gets a correction proportional to a;*^™'^^-'", where x 
is the common value of all |xfc|. Since the bare m + 1-particle cumulant Q^! scales like 2;-(™+i)(i+") (cf. Sec. ^), the 
anomalous exponent is again exactly canceled by the renormalization. This is important, since one needs a to scale 
to zero not only in the self-energy, but also in the interaction vertices, in order for the low-energy fixed point to be 
asymptotically free. 
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FIG. 13. d \og[Fm{m l,V)\/d I vs I for tti = 4 (solid line) and m = 2 (dashed), with a logarithmic Fade approximation 
d \og[To(l)]/ dl = P„il)/Qkil) with n = fc = 19. 
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